Small Angle Approximation for non parallel plate 
capacitors with Applications in Experimental 
Gravitation 

B R Patla 

Smithsonian Astrophysical Observatory, 
Harvard-Smithsonian Center for Astrophysics, 
60 Garden St, Cambridge, MA 02138, USA 

E-mail: bpatla@cfa.harvard.edu 

Abstract. An approximate analytical formula for the capacitance of a non-parallel 
plate capacitor with small values of inclination angles and distance separations of the 
plates is presented. Most applications involving position sensing that are employed 
in precision experiments, micromechanical devices, etc often use a parallel plate 
approximation for modeling the capacitance. The analytical approximation presented 
here is based on the more general, but hard to implement formalism for small angles 
given in Xiang (2005), and may be used directly to compute the forces acting on the 
proof masses by taking the derivative of the capacitance. 

Effects of varying the inclination angle, plate dimension and separation on the value 
of the capacitance per unit length are analyzed in detail. We use the formula derived in 
this paper to compute the acceleration of the proof mass due to the presence of sensing 
electrodes (used for positioning the proof masses) with an assumed tilt of ~ 10 -6 rad as 
applicable to SR-POEM: an experiment aiming to test the weak equivalence principle 
- a fundamental postulate of general relativity [7J [5] . 



PACS numbers: 04.80.Cc, 07.05.Fb, 07.05.Tp 
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1. Introduction 

Capacitive sensing is one of the most commonly used techniques in science and 
engineering for position measurements ranging from a few millimeters up to slightly 
less than a nanometer. In experimental physics, often a parallel plate approximation 
is used for positioning the proof masses of gravitational wave detectors [HI H], precision 
experiments like SR-POEM that intends to test the weak equivalence principle (WEP), 
torsion balance experiments [8], etc. 

More recently, the effects of assuming a non-parallel plate capacitor model for the 
switches used in micro electro mechanical systems (MEEMS) [TT] and differential micro 
accelerometers used in micromechanical machining [9] have both concluded that using 
the simplified parallel plate model underestimates the capacitance by a few percent. 

In the case of WEP test, before the actual measurement of the precise positions 
of the proof masses in free-fall commences, the proof masses have to be positioned 
accurately using electrostatic force within the housing. If the proof masses are let non- 
parallel with respect to the sensing and forcing electrodes after positioning, the resulting 
configuration will introduce a component of force (on the proof masses) along the 
direction of the WEP measurement which will be slightly more than the corresponding 
configuration mimicking a parallel plate configuration. A rough estimate of the force can 
be obtained by taking the derivative for the capacitance of two parallel plates, eoA/d, 
where A is the area of the plate, d the distance separating the plates and eo is the 
vacuum permittivity, by neglecting the relative tilt of the plates. 

The general formula for non-parallel plates as given in [12] involves elliptic integrals 
that have to be evaluated numerically. For tiny distance separations and small angles 
that are of interest in experimental gravitation, the integrals evaluate to infinity even 
with reasonable numerical precision because of the very nature of elliptic integrals 
being singular at unity and thereby requiring one to impose cut-offs, resulting in an 
over determination of the capacitance. In order to obtain reasonable estimates of the 
capacitances to first order approximation, a convenient formula for capacitance that 
establishes a functional relationship between angle and distance separating the plates 
is more useful. Moreover, the force between the plates is directly proportional to the 
derivative of the capacitance and the closed form approximate solution of the capacitance 
is useful for computing the force directly and more accurately. 

In section [2] we provide a brief and transparent overview of the results in [12] , whose 
results we will refer to as semi-analytical from hereafter throughout this paper because of 
the elliptic integrals involved in the final form that have to be evaluated numerically. To 
make clear distinction, the results derived in section [3] will be dubbed as the small angle 
approximation, which is valid only for small values of inclination angles and separations 
(much smaller than the linear dimensions of the plates). The result of this section is 
compared with finite element analysis (FEA) model and also the semi-analytic formula, 
section [4] delineates the effects of the relative orientation of the plates using the semi- 
analytic formula, because there is no difference in the value of capacitances to first order 
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using the formula for small angle approximation. In section [5] we will use our small angle 
approximation formula to study the spurious WEP (acceleration) signal in SR-POEM , 
followed by a brief discussion in section [6j 

2. Capacitance of two non-parallel plates 

We briefly introduce the formalism developed in [12] to set the stage for arriving at an 
approximate solution for the capacitance per unit length in the next section. Consider 
two infinite strips extending through the z direction and subtending an angle <fi between 
them, see figure flj The coordinates of the center of the upper plate are (Ax, Ay) from 




Figure 1. Plates of the non parallel plate capacitor have widths Liand L 2 respectively 
with their lengths extending infinitely along the z direction. The geometric centers of 
the plates coincide along the horizontal axis (Ax = 0, need not be true in general) 
and are separated by a distance Ay = h along the y axis. The extensions of the plates 
that intersect at a point O are the lengths r% and r 2 . The plates make an angular 
separation <p (measured counter clockwise from the x axis). Points A, B, C and D are 
the end points of the plates in the xy plane. 



the center of the bottom plate. The configuration is subjected to successive multiple 
coordinate transformations as given in Table [TJ The idea is to find a coordinate system 
in which the non-parallel plates of unequal widths are parallel to each other and also 
have equal widths. The capacitance of such a configuration is to be identified with that 
of an ordinary parallel plate capacitor, albeit the linear dimensions corresponding to the 
aforesaid configuration being complex functions. 

The initial configuration in the old complex z plane is transformed on to the new 
t plane such that the ends of the plates are fixed at +1 and —1 respectively. These 
constraints determine the functional form of the constants to be 

Cl ~ (Li + nyi* + (l 2 + r 2 yi* and 

(Lx + riff - (L 2 + r 2 fl* , 
°° {L 1 + r i y/'f> + {L 2 + r 2 y/^ U 
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Table 1. Set of successive coordinate transformations of the non-parallel plates 
Coordinates^ 

Plane A B CD Transformation 

z ne* (ri+Li)e* r 2 r 2 + L 2 

t (a, 0) (-1,0) (/?,0) (1,0) t = c l z v ^ + c a *, 0^tt 

C (-00,0) (0,0) (i/fc 2 ,o) (1,0) C - 

u iiT' K + iK' K Schwarz-Christoffel 

t End points of the cross-section of the upper and the lower plates, figure 1. 

* Constants c\ and cq are determined by imposing the constraint — x coordinates of 

B and D be — 1 and 1 respectively. 



where L% and L 2 are the widths of the plates respectively. The transformation into the 
C plane (see Table [T]) yields coordinate distances that are appropriate for the Schwarz- 
Christoffel transformation [10J that is to be invoked next, which maps the interior of 
the plates on to the area of a rectangle. Solving the equation (k 2 — 1=0, yields 

where 

-2rl'* + (Li + n)^ - (L 2 + rgff 

(Li + n)*/* + (L 2 + r 2 Y^ ^ 

2r ^ + (L 1 + r 1 )^-(L 2 + r 2 )^ 

^ (L 1+ri )*/* + (L 2 + 7 . 2 )ir/* • 1 J 

The Shwarz-Christoffel transformation maps the upper half of the complex plane into 
the interior of a rectangle. 

/ \ — — 1 



^ ^(c-o^cc-i)^- 1 ^-^) (4) 



u = A 1 k dC + A . (5) 

io VC(i-0(i-Cfc 2 ) 

Substituting ( = p 2 and imposing the constraints at points -B(A) = 0) and = 1) 

respectively, we identify 

u - s K(k] = 1 7w = m = m dp and (6) 

where is the complete elliptic integral of the first kind of modulus with 

K(k') = K'(k) and k' 2 = 1 - k 2 . From equation g, 

/ (l + 0X1-/3) 

" = V(i-«)(i + /3)- (8) 
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The capacitance per unit length, using the parallel plate formula, is 

where k- m and k out represent the values of k (equation ([2])) for the interior and exterior 
of the plates corresponding to the angles and 2n — <p in equation 

3. Approximate formula for small angles 



Values of the elliptic integrals in equation (|9j) have to be evaluated numerically. For 
very small values of <fi and h/L (see figure [T|, even after allowing a reasonably high 
numerical precision, we encounter indeterminate values for the integrals and so one has 
to impose cut-offs very close to unity. These numerical manipulations often introduce 
errors that result in over or under determination of the capacitances based on the 
numerical algorithms used. In order to compute the force as a function of distance 
separating the plates and the included angle between them, we have to numerically 
evaluate the value of capacitances for a range of angles and separations followed by 
numerical differentiation. We will use a simplified model for the capacitor plates by 
demanding their lengths be equal, L\ « L 2 ~ L. From A AOC, 

n « -— - , r 2 w — (10) 

sin <p sin <p 

T h h 2 s'm(5 2 - <p) . . 

r l + L l ^^— , r 2 + L 2 ~ — 11 

sin (p sm (p 

h 1 = h-^<l>, h 2 = h+^<P (12) 
From equations ((8| and p 



where 5\ = ir/2 — (p,5 2 = n/2 + <p and L\ = L 2 = L. Using equations (10) and (11) in 



equation (13) yields 

k't = tanh tanh (tt0 + . (14) 

equation (14) is valid when h/L > (ft and <C 1, necessary conditions for small 
angle approximation. Note that, the usual condition for parallel plate capacitor, 
L/h ^> 1(~ 10 is sufficient in most cases), must also be satisfied independently. Letting 

— = a and ^+^- = 6 (15) 



where a,&> 1, yields 

v^-l-^— (1 + e- 2 ^). (16) 



e -2a 
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The relationship between the elliptic integrals and their complementary modulus may 
be expressed in terms of the nome 
K'(k) 1 



logg, 



(17) 



K(k) 7T 

where q is the nome which may also be expressed as a functional of theta functions 
denoted by A (see appendix for a derivation). 

1 - \/k 



A ; 



e~ 2a (1 - 7T0) (4 - 2e~ 2a (1 



7T 



18) 



'"" 2(1 + Vk>, 

ignoring the second and higher order terms in the expansion of the exponential term. 

In order to evaluate the exterior capacitance for the configuration in figure [l] by 
using equation (JoJ) , we need an expression for k! outside. Replacing <p — > 27r — cf> in 



equation (jX3h one obtains k' 2 OVLt = k'^A^^w-^ &'<mt ~ Lrft/Ah, yielding 



A 



out 




(19) 



where we replaced k' in with k' out in equation (pL8|) . The capacitance per unit length, to 
first order approximation in <p is 

C (i) = _ eo 



7i 



(log A in + log A 



out) 



(20) 



Including higher order terms in A and using equations (18) and (19) in equation (20) 
we obtain an expression for the capacitance 



C 



h 



71 



log (1 - 7T0) + log 



+ 



eo 



7T 

£o 

71 



log 4 ^2 - (1 - 7T0) exp 
-^[log(l + 2Af n + ---)+log(l + 2A 





+ ... 



(21) 



The two terms in the first line of equation (21 ) comprises the equivalent of the formula 



for the parallel plate capacitor (in the limit of <p — > 0) and the last two terms represent 
the first order fringe capacitance and the rest of the higher order terms. The result of 
this section is compared with numerical simulations (FEA) and [T2j , and are plotted in 
figure [2j For the FEA, when = 0, the value of the capacitance match with the standard 
parallel plate capacitor formula to within 2%. As the tilt of the upper plate is increased, 
the error in the estimated value of the capacitance increases due to the leaking of the 
fringe capacitance into the bounding box. The small angle approximation matches with 
the FEA more than it does for the semi-analytical formula due to over estimation of 
the capacitance introduced as a result of cut-offs imposed on the numerical estimation 
of the elliptic integrals. Note that equation ^ contains the ratio of elliptic integrals 
and that is precisely the reason why the value of the capacitance is over determined 
instead of being otherwise. However, for values of ~ 10~ 2 and above, the small 
angle approximation is inaccurate and the semi-analytic approach is more reliable as 
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the elliptic integrals now yield finite values. Therefore, the approximate solution for 
small angles is a rather convenient way to estimate the capacitance and the force for 
very small angles that we often encounter in experimental gravitation. The capacitance 
is now a function of the distance separating the plates and also the inclination angle. 



4. Capacitance as a function of rotation angle: Two simple cases 

Although for small angles and separations the choice of the rotation axis of the upper 
plate (figure [T]) has no effect on the capacitance to first order in 0, it matters for large 
values of inclination angles. This is due to the difference in the values of the effective 
separation between the plates corresponding to the choice of rotation axes. The semi- 
analytic approach outlined in [12] is used to investigate two cases corresponding to two 
different axes of rotation, for inclination angles that are greater than 10~ 2 rad. 

We consider two simple cases wherein the axes of rotations of the upper plate are 
1) a line passing through the middle (along the z direction in figure [T]) of the upper 
plate and 2) an edge of the upper plate directly above a corresponding edge of the 
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Figure 2. Capacitance corresponding to a pair of plates each with sides 1 mm 
separated by a distance of 0.1 mm for angles 10~ 4 , 10~ 3 , 10~ 2 and 0.1 rad. The data 
points are obtained by using three different methods: 1) Semi-analytical from [TJ] 2) 
Finite Element Analysis (FEA) 3) Small angle approximation. The upper plate is 
tilted about its center with the bottom plate fixed. The curves connecting the data 
sets are quadratic best fits corresponding to the three cases. The error in the FEA is 
less than 5% for the maximum values of <f> in the figure. Fringing capacitance is not 
included in this plot. 
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lower plate. We will study the variation of capacitance with angle, plate separation, 
and relative dimensions of the plates. 

4-1. Rotation about the centerline of the upper plate 

We introduce a coordinate system described in figure [T] with 



sin . 




ri = ^7 " ^ ; r lmin = J( I ^ + Ax) + Ay 2 - ^, (22) 



r 2 = 7—^7 - Ax - ^ ; max = sin 1 ( Ay ) , (23) 

tan0 2 y^ +rimia J 

where r lmin is the minimum distance of the nearest edge of the upper plate from the 
point of intersection of the plates and max is the maximum value allowed for the rotation 
angle, with all other notations retaining the same definitions as that of in the caption 
of figure [l] in section |2} Using equations (22) and (23) in equation ^ we proceed to 



study the capacitance as a function of included angle and relative linear dimensions of 
the plates for distance separations of 1 mm, 5 mm and 1 cm respectively. The chosen 
separations are the typical values encountered in SR-POEM - the spacing between the 
proof mass and housing - discussed at length in section [5] 

For plates of equal widths, the effective separation decreases for the left half of 
the configuration and thereby increasing the capacitance. The total capacitance is the 
sum of the capacitances as if they are in parallel yielding a smaller effective separation 
(product of the distance separations in the formula for two capacitances in parallel) 
resulting in a higher capacitance as the angle is increased. As the offset along the 
x direction (Ax) is increased, the plates overlap less if their linear dimensions are 
held constant, thereby reducing the capacitance. As the included angle is increased 
the effective separation of the left half increases, contributing even more to the total 
capacitance. These results are plotted in figure [3] 

4-2. Rotation about an edge of the upper plate 

We introduce a coordinate system by redefining 

n = - — - and r 2 = -. (24) 

sin tan 

where Ay = h. The plate separation is now the distance separating the aligned edges 

of the two plates of equal widths and being the inclination of the upper plate with 

respect to the horizontal lower plate. The variation of capacitance as a function of and 

relative linear dimensions are shown in figure |4j The capacitance decreases with as 

opposed to the case when the axis of rotation was directly above the center of the lower 



plate( section 4.2). This is because the effective distance is now directly proportional to 
0. But the functional dependence of capacitance on the linear dimensions (L\ and L2) of 
the plates exhibit a similar pattern to that of the configuration described in section |4.2[ 
see figure |4j 
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Figure 3. Variation of capacitance per unit length for a pair of plates, each with width 
5 cm for distance separations corresponding to 1 mm, 5 mm and 1 cm respectively, (a) 
As a function of included angle <j>. (b) As a function of plate offset Ax. The included 
angle is fixed at 1CT 2 rad. 




0.05 0.10 0.15 0.20 1 2 3 4 5 

</> (rad) i2 ( cm ) 



Figure 4. Variation of capacitance per unit length for distance separations 
corresponding to 1 mm, 5 mm and 1 cm respectively. The dashed line includes the 
fringing capacitance, which is less than a percent, (a) The widths of the plates are 
fixed at 5 cm each. (b)The upper plate is set at an angle of 10 -2 rad with respect to 
the lower plate and L\ = 5 cm. 
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5. Test of the weak equivalence principle: A case study 

Testing of the equivalence principle as in SR-POEM amounts to the determination of rj 
parameter (measure of relative acceleration of two test masses of different materials) to 
an accuracy of 10~ 17 m/s 2 [7\. A laser gauge monitors the positions of the test masses, 
as the proof masses inside the vacuum chamber become part of a high finesse Fabry- 
Perot cavity Before the measurements are made the test masses are properly aligned 
electrostatically by applying a drive voltage. When the actual distance measurements 
are performed, the capacitive sensing is turned off. However, the surface potential 
difference between the test mass and the sense electrodes that are used to align the test 
masses will generate a force along the measurement direction and produces a spurious 
acceleration signal. Many factors are responsible for generating a surface potential 
difference; patch effect being the most important one and that could be as high as a few 
milli volts [31 12]. 

Even if the proof masses are properly positioned before the measurement, i.e. to say 
the sides of the proof masses are exactly parallel to the walls of the housing, there will 
still be a force acting on the proof mass in the direction perpendicular to the direction 
of the position measurement, and hence is of little concern to us. This is because a side 
of the proof mass together with a side of the housing will constitute the equivalent of a 
parallel plate capacitor. However, imperfections in machining or the rotational inertia 
(as soon as the drive voltage is turned off) of the proof mass will cause a small tilt of 
the proof mass with respect to the housing walls. We will assume a maximum value of 
10 -6 rad for this misalignment (tilt). A small component of the force due to this tilt 
will result as an acceleration along the WEP measurement axis. The force acting on 
the test mass is 

F = -±V 2 VC(<j>,h) (25) 

where C ((/>, h) is the capacitance and V is the (surface) potential difference across the 
plates. The resultant force acting on the test mass, for a relative inclination angle of 
with any given side of the housing, is 



7T 



,V 2 sm<t> (26) 
1- 7i(f) n^{s 2 -l)J ^ V ; 

where s = (//4/i) 1//2 , / being the linear dimensions of the test masses and h is the distance 
separating the test mass and the sensing electrode. If there is a sensing electrode on the 
bottom (although it is not a good idea to have one), the force is 



2s 3 ( 4s - . y J . ) V 2 + ( - , ) V'cos^. (27) 

1 7r{s 2 -l)J V 1 "^ ttv^s 2 -!),' 



In deriving the above results we have neglected the fringing capacitance in equation (21 ). 
The variation of the capacitance (with and without fringe effects) using small angle 
approximation derived in section [3] is compared with the parallel plate case in figure |5j 
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Figure 5. Variation of capacitance per unit length as a function of distance. Cf is the 
fringing capacitance, Cp corresponds to the well-known formula for the parallel plate 
case, €qL/1i. Assumed plate dimensions are 5 cm with an inclination angle of 10 -6 rad 
and a potential difference of 10~ 3 V between the plates. 

For SR-POEM, estimate the force on the test masses surrounded by sensing 
electrodes of sides 5 cm and separated by a nominal distance 5 mm from the test mass, 
allowing for a tilt of = 10~ 6 rad and a potential difference of a milli volt. 

(i) A plate facing a side of the test mass with the measurement direction perpendicular 
to the sides of the test mass will produce a force 4.48 x 10 -21 N 

(ii) A bottom plate of similar dimensions yields a higher force 4.64 x 1CT 15 N. This 
configuration is usually avoided. 

(iii) A parallel plate capacitor formula yields a force 8.85 x 10~ 16 N 
6. Discussion 

We have studied the effects of non-parallel plate capacitors in the context of experimental 
gravitation where the angles are very small. An approximate formula to compute 
the capacitance in this regime is thoroughly analyzed along with the semi-analytic 
approach given in [12J and finite element methods. For very small angles the small 
angle approximation matches well with the finite element simulations (with maximum 
2% error in the appropriate angular regime). Around and after cf) ~ 10~ 2 rad the 
semi-analytical formula is a closer match. 

In conclusion, we have shown that the small angle approximation for estimating the 
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capacitance of non-parallel plate capacitors is a very good alternative to computationally 
expensive numerical methods and semi-analytical approaches. For applications in 
experimental gravitation it is easier to compute the forces, acceleration and also estimate 
other noise sources due to positioning of test masses using our approximate solution. 
Our analysis successfully demonstrates that for an assumed surface potential difference 
of 1(T 3 V or less SR-POEM will easily meet it's mission goal of 10 17 m/s 2 . 

Appendix A. Appendix 

The nome q is defined as 

q = exp(-nK'/K). (A.l) 
The elliptic functions may be represented in terms of the theta functions 

oo 

^(z,q) = Y,(-lT~k (nH)2 e {2n+1)lz 

n=— oo 

oo 

<& 2 (z,q) = q {n+1 i )2 e {2n+1)lz 

n=— oo 

oo 

<& 3 (z,q) = q n2 e 2niz (A.2) 



n=— oo 
oo 

2niz 



n=— oo 

■&j(z + 7r, q) = =F "&j(z, q) negative sign for j = 1, 2 
$j(z + 7T7, q) = =p N$j(z, q) negative sign for j = 1, 4 

where N = q~ l e~ 2lz . A convenient way to obtain the ratio K' /K, is by taking the 
logarithm of the nome q. Using the theta functions above, it is possible to arrive at a 
differential equation, the solution of which is the elliptic integral pU]. The form of the 
modulus is 

k = U^q) 2 /M0,q) 2 ^ 2 2 /el (A.3) 
where we have adopted a new notation for the theta function corresponding to z — 0. 



Below we write down the series expansions using equation(A.2) for z = 

9 2 = 2q 1 ^ + 2g 9 / 4 + 2g 25 / 4 + 2g 49 / 4 + 2g 81 / 4 + 2g 121 / 4 + • • • 

9 3 = 1 + 2q + 2q 4 + 2q 9 + 2q w + 2q 25 + 2q m + ■■■ (A.4) 
e 4 = 1 - 2q + 2q 4 - 2q 9 + 2q m - 2g 25 + 2g 36 + • • • 



noting that Q\ = 0. Also, since Q\ + Q\ = 9 A , using equation (A.3) we have k' = Q 2 JQ\- 
The set of equations (A.4) may be combined to obtain a power series in q. 
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q + q 9 + q 25 + q i9 + --- .. r . 

(A. 5) 



1 + 2g 4 + 2g 16 + 2g 25 + 2g 36 + 



Inverting the series in equation (A.5), we obtain 



q = A + 2A 5 + 15A 9 + 150A 13 + 1725A 17 + • • • (A.6) 
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